A Local Model for Angular Momentum Transport in Accretion Disks 
Driven by the Magnetorotational Instability 
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We develop a local model for the exponential growth and saturation of the Reynolds and Maxwell 
stresses in turbulent flows driven by the magnetorotational instability. We first derive equations 
that describe the effects of the instability on the growth and pumping of the stresses. We highlight 
the relevance of a new type of correlations that couples the dynamical evolution of the Reynolds and 
Maxwell stresses and plays a key role in developing and sustaining the magnetorotational turbulence. 
We then supplement these equations with a phenomenological description of the triple correlations 
that lead to a saturated turbulent state. We show that the steady-state limit of the model describes 
successfully the correlations among stresses found in numerical simulations of shearing boxes. 
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Since the early days of accretion disk theory, it has 
been recognized that molecular viscosity cannot account 
for a number of observational properties of accreting ob- 
jects. Shakura and Sunyaev [l| introduced a parametriza- 
tion of the shear stress that has been widely used since. 
Much of the success of their model lies on that many 
disk observables are determined mostly by energy bal- 
ance and depend weakly on the adopted prescription . 
However, this parametrization leaves unanswered funda- 
mental questions on the origin of the anomalous trans- 
port and its detailed characteristics. 

Strong support for the relevance of magnetic fields in 
accretion disks arose with the realization that differen- 
tially rotating flows with radially decreasing angular ve- 
locities are unstable when threaded by weak magnetic 
fields 0| ■ Since the discovery of this magnetorotational 
instability (MRI), a variety of local jj, |a] and global 
0- ll| numerical simulations have confirmed that its 
long-term evolution gives rise to a sustained turbulent 
state and outward angular momentum transport. How- 
ever, global simulations also demonstrate that angular 
momentum transfer in turbulent accretion disks cannot 
be adequately described by the Shakura & Sunyaev pre- 
scription, fn particular, there is evidence that the tur- 
bulent stresses are not proportional to the local shear |i| 
and are not even determined locally [l(ij|. 

Some attempts to eliminate these shortcomings have 
been made within the formalism of mean-field magneto- 
hydrodynamics [Tl| . This has been a fruitful approach 
for modeling the growth of mean magnetic fields in differ- 
entially rotating media 12]. However, it has proven diffi- 
cult to use a dynamo model to describe the transport of 
angular momentum in these systems This is es- 

pecially true in turbulent flows driven by the MRI, where 
the fluctuations in the magnetic energy are larger than 
the mean magnetic energy, and the turbulent velocity 
and magnetic fields evolve simultaneously Q. 

In order to overcome this difficulty various approaches 
have been taken in different physical setups. Blackman 



and Field [lfij derived a dynamical model for the nonlin- 
ear saturation of helical turbulence based on a damping 
closure for the electromotive force in the absence of shear. 
On the other hand, Kato and Yoshizawa [lH Il7j and 
Ogilvie derived a set of closed dynamical equations 
that describe the growth and saturation of the Reynolds 
and Maxwell stresses in shearing flows in the absence of 
mean magnetic fields. 

In this Letter, we relax some of the key assumptions 
made in previous works and develop the first local model 
for the dynamical evolution of the Reynolds and Maxwell 
tensors in turbulent magnetized accretion disks that in- 
corporates explicitly the MRI as their source. 

The equation describing the dynamical evolution of the 
mean angular momentum density, I, of a fluid element in 
an accretion disk with tangled magnetic fields is 



d t l + V-(Zt>) = -V-(rJF) 



(1) 



Here, the over-bars denote properly averaged values, v 
is the mean flow velocity, and the vector F character- 
izes the flux of angular momentum. Its radial compo- 
nent, T r = R rt j} — M rt p, is equal to the total stress acting 
on a fluid element. The stresses due to correlations in 
the velocity field, R r ^, = (pSv r Sv^), and magnetic field 
fluctuations, M ri j> = (SB r SB^,} /4tt, are the Reynolds and 
Maxwell stresses, respectively. A self-consistent accretion 
disk model based on the solution of the equations for the 
mean velocities requires a closed system of equations for 
the temporal evolution of these mean stresses. 

In a recent paper j20j . we identified the signature of 
the magnetorotational instability in the mean Maxwell 
and Reynolds stresses due to correlated fluctuations 
of the form Sv = [5v r (t, z), Sv^(t, z), 0] and SB — 
[SB r (t, z), SB<p(t, z),0] in an incompressible, cylindrical, 
differentially rotating flow, threaded by a mean verti- 
cal magnetic field, B z . We demonstrated that a num- 
ber of properties of the mean stresses during the initial 
phase of exponential growth of the MRI are approxi- 
mately preserved in the saturated state reached in local 
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three-dimensional numerical simulations. As a first step 
in our calculation, we aim to derive dynamical equations 
for the mean stresses that describe this result. 

In order to work with dimensionless variables we con- 
sider the characteristic time- and length-scales set by 
1/fio an d ^Az/^o- Here, Q and v\ z — B z /y/4npo denote 
the local values of the angular frequency and the (verti- 
cal) Alfven speed in the disk with local density po at the 
fiducial radius ro . The equations governing the local dy- 
namics of the (dimensionless) MRI-driven fluctuations in 
Fourier space, are then given by [20| 



d t 6v r = 28v < t > + ik n Sb r , 

dtSvff, = (q - 2)Sv r + i/c„(5&0 . 

dtSb r = ik n 8v r , 

dtSbcf, = -q5b r + ikndvcf, , 



(2) 
(3) 
(4) 
(5) 



where Svi and Sbi stand for the Fourier transform 
of the dimensionless physical fluctuations 5vi/vAz and 
8Bi/B z . The wavenumber, k n , denotes the mode with 
n nodes in the vertical direction and the parameter 
q = — dhiCl/dhir\ ro is a measure of the local shear. 

Using the fact that the modes with vertical wavevec- 
tors dominate the fast growth driven by the MRI, we 
obtain a set of equations to describe the initial exponen- 
tial growth of the Reynolds and Maxwell stresses. We 
start from the equations for the fluctuations @~© and 
use the fact that the mean value of the product of two 
functions / and g, with zero means, is given by |20| 



(fg)(t) = 2j2 Re[f(k n ,t)g*(k n ,t)] 



(6) 



71=1 



By combining different moments of equations we 
obtain the dimensionless set 

dtRrr = 4R r<f) + 2W r4> , (7) 

d t R r< p = {q-2)R rr + 2R H -W rr + W H , (8) 

dtR^ = 2(g-2)i? r0 -2WV, (9) 

d t M rr = -2W r</) , (10) 

d t Mr<t> = -qMrr + Wrr-W^, (11) 

dtM^ = -2qM r4> + 2W^ r , (12) 
where we have defined the tensor Wij with components 



W rr = (Sv r Sj r ) = (Sb^SuJff,) 

W r<i> = (SVrSjtj,) = -{SbrSUtf,} 

W^r = (5v<f,5j r ) = -{Sb^SUr) 

= (dv^Sjj,) = (Sb r 5u)r) 



(13) 
(14) 
(15) 
(16) 



Here, 5ji and StOi, for i = r,cj>, stand for the compo- 
nents of the induced current Sj = VX(56 and vorticity 
5u) = Vx<5i> fluctuations. Note that the components of 
the tensor Wij are defined in terms of the correlations 



between the velocity and current fields. However, for the 
case under consideration, these are identical to the cor- 
responding correlations between the magnetic field and 
vorticity fluctuations (eqs. |T3]-[ro]). 

Equations l(7|l- ljT2|) show that the MRI-driven growth 
of the Reynolds and Maxwell tensors can be described 
formally only via the correlations Wij that connect the 
equations for their temporal evolution. Note that, in 
contrast to the correlations (SviSuik) and (SbiSjk), that 
appear naturally in helical dynamo modeling and trans- 
form as tensor densities (i.e., as the product of a vec- 
tor and an axial vector), the correlations (Svidjk) and 
(SbiSuik) transform as tensors. Moreover, the tensor Wij 
cannot be recast in terms of the cross-helicity tensor 
Hij = (SviSbj) because, for the unstable MRI modes, 
the ratio Hij /Wij approaches zero at late times. Neither 
can Wij be expressed in terms of the turbulent electro- 
motive force (5vx8b), as the latter vanishes under our 
set of assumptions, implying that no mean magnetic field 
is generated. 

In order to describe the MRI-driven exponential 
growth of the stresses in equations (|T|)- ((T2")l we need to 
write an additional set of dynamical equations for the 
evolution of the tensor Wij . Using appropriate combina- 
tions of different moments of equations ©-I© we obtain 

OtWrr = qW r <t> + 2W+r + {l^f R r4> - (k™) 2 M rtj} , (17) 

d t W r4l = 2W U - (k? r ) 2 R rr + {k^) 2 M rr , (18) 
dtW^r = (q- 2)W rr + qWn + Q$+?Rm> ~ Q&fMn, , (19) 
dtW^ = {q- 2)W r4> - (k^) 2 R r(j> + {k^fM r4> , (20) 

where we have defined the set of mean wavenumbers 



J2 k 2 n M5bi5b*] 



R\2 — n=l 



Re[SviSvj 



(k 



M\2 _ n=l 



Re[5biSbj 



The system of equations GJ-lO and |TT|— (JSOJ de- 
scribes the temporal evolution of the stresses during the 
exponential growth of the MRI in a way that is formally 
correct, with no approximations. Motivated by the simi- 
larity between the ratios of the stresses during the expo- 
nential growth of the MRI and during the saturated tur- 
bulent state [2(J, we propose to use the right-hand sides 
of equations f?j)- fT2"|) and fT7jl - (|2l1|l as a local model for 
the source of turbulence in MRI-driven magnetohydro- 
dynamic flows. Of course, in the turbulent regime, the 
various average wavenumbers k^' will depend on the 
spectrum of velocity and magnetic field fluctuations. As 
the lowest order model for these wavenumbers we choose 



i, j = r,<t>, (21) 
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FIG. 1: Correlations between the Maxwell and Reynolds stresses 
and mean magnetic energy density at saturation in MRI-driven 
turbulent shearing boxes Q. The lines show the result obtained 
with our model in the asymptotic limit for £ = 0.3. 



FIG. 2: The mean magnetic energy density in terms of a sat- 
uration predictor found at late times in numerical simulations of 
MRI-driven turbulence Q . The line shows the result obtained with 
our model in the asymptotic limit for f = 0.3 and £ = 11.3. 



where £ is a parameter of order unity and k max corre- 
sponds to the wavenumber at which the growth rate of 
the fluctuations reaches its maximum value, 7 max = <?/2. 

By construction, the set of equations l(7)l- i|12|) and 
(I17|l l|20|) leads to the expected exponential growth of 
the mean stresses driven by modes with wavevectors per- 
pendicular to the disk midplane. These modes, however, 
are known to be subject to parasitic instabilities, which 
transfer e nerg y to modes in the perpendicular (k r ,k^,) 
directions |l9j . The initial fast growth experienced by 
the stresses will eventually be slowed down by the com- 
bined effects of non-linear couplings between modes and 
of dissipation at the smallest scales of interest. 

The terms accounting for these interactions would ap- 
pear in the equations for the stresses as triple correla- 
tions between components of the velocity and magnetic 
fields, i.e., they would be of the form (SviSvjSbk) and 
(SvidbjSbk) ■ These types of non-linear terms have also 
been considered by Ogilvie |18| . who proposed scalings 
of the form 

{SviSvjSbk) ~ {SviSvj) (SbM 1 / 2 ~ M^r.. . (22) 

In the absence of a detailed model for these correla- 
tions and motivated again by the similarity of the stress 
properties during the exponential growth of the MRI and 
the saturated turbulent state [2(|, we introduce a phe- 
nomenological description of the non-linear effects on the 
evolution of the various stresses. In particular, denoting 
by Xij the ij-component of any one of the three tensors, 
i.e., Rij, Mij, or Wij, we add to the equation for the 
temporal evolution of that component the sink term 



dt 



sink 




(23) 



Here, M/2 = (M rr + M<^)/2 is the mean magnetic en- 
ergy density in the fluctuations and Mq is a parameter. 



Adding the sink terms to the system of equations JJJ 
B12J1 and Ijl7|l -I|20 |l leads to a saturation of the stresses af- 
ter a few characteristic timescales, preserving the ratio of 
the various stresses to the value determined by the expo- 
nential growth due to the MRI and characterized by the 
parameter £. These non-linear terms dictate only the sat- 
urated level of the mean magnetic energy density accord- 
ing to lim^oo M = r 2 Mo, where T is the growth rate for 
the stresses, which, in the case of a Keplerian disk, with 

q = 3/2, is given by T 2 = 2 Ul + 15C 2 - (l + 15C 2 /8) . 

We infer the dependence of the energy density scale 
Mo/2 on the four characteristic scales in the problem 
f^Oj H (the vertical length of the box), po, and vaz using 
dimensional analysis. We obtain Mq/2 oc poH 5 QqV a ~ 5 
which leads, with the natural choice S = 1, to Mq/2 = 
^PqHQqvaz, where we show the dimensional quantities 
explicitly and introduce the parameter £. 

Our expression for Mo describes the same scaling be- 
tween the magnetic energy density during saturation and 
the various parameters characterizing the disk found in 
a series of shearing box simulations threaded by a finite 
vertical magnetic field and with a Keplerian shearing pro- 
file 0. By performing a numerical study of the late-time 
solutions of the proposed model, we found a unique set 
of values (C,£) such that its asymptotic limit describes 
the correlations found in these numerical simulations. 

Figure^shows the correlations between the ^-compo- 
nents of the Maxwell and Reynolds stresses and the mean 
magnetic energy density found during the saturated state 
in numerical simulations 0. In our model, this ratio de- 
pends only on the parameter £ and the shear q (held fixed 
at q = 3/2 in the simulations). It is evident that, in the 
numerical simulations, the ratios of the stresses to the 
magnetic energy density are also practically independent 
of any of the initial parameters in the problem that de- 
termine the magnetic energy density during saturation 
(i.e., the x-axis in the plot). This is indeed why we re- 
quired for our model of saturation to preserve the stress 
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ratios that are determined during the exponential phase 
of the MRI. Assigning the same fractional uncertainty 
to all the numerical values for the stresses, our model 
describes both correlations simultaneously for £ = 0.3. 

Figure [21 shows the mean magnetic energy density in 
terms of a saturation predictor found in numerical simu- 
lations Q. For £ = 0.3 and assigning the same fractional 
uncertainty to all the numerical values for the stresses, 
we obtain the best fit for £ = 11.3. These values for the 
parameters complete the description of our model. 

In summary, in this Letter we developed a local model 
for the evolution and saturation of the Reynolds and 
Maxwell stresses in MRI-driven turbulent flows. The 
model is formally complete when describing the ini- 
tial exponential growth and pumping of the MRI-driven 
stresses and, thus, satisfies, by construction, all the math- 
ematical requirements described by Ogilvie [lij . Al- 
though it is based on the absolute minimum physics 
(shear, uniform B z , 2Z?-fluctuations) for the MRI to be at 
work, the model is able, in its asymptotic limit, to recover 
successfully the correlations found in three-dimensional 
local numerical simulations 

Finally, the local model described in this Letter con- 
tains an unexpected feature. The mean magnetic field in 
the vertical direction couples the Reynolds and Maxwell 
stresses via the correlations between the fluctuations in 
the velocity field and the fluctuating currents generated 
by the perturbations in the magnetic field. These second 
order correlations, that we denoted by the tensor Wij, 
play a crucial role in driving the exponential growth of 
the mean stresses and energy densities observed in nu- 
merical simulations. To our knowledge this is the first 
time that their relevance has been pointed out in either 
the context of dynamo theory or MRI-driven turbulence. 
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